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ABSTRACT 

Using the MOPED algorithm we determine non-parametrically the Stellar Mass Function 
of 96,545 galaxies from the Sloan Digital Sky Survey data release one. By using the re- 
constructed spectrum due to starlight we can eliminate contamination from either emission 
lines or AGN components. Our results give excellent agreement with previous works, but 
extend their range by more than two decades in mass to 10^'^ < Ms/h~'^MQ < 10^^. 
We present both a standard Schechter fit and a fit modified to include an extra, high-mass 
contribution, possibly from cluster cD galaxies. The Schechter fit parameters are (/>* ~ 
(7.7 ±0.8) X IQ-^h^MpC^^, M* = (7.53 ±0.04) x W^^H-^Mq anda = -1.167 ±0.004. 
Our sample also yields an estimate for the contribution from baryons in stars to the critical 
density of r^h* /i = (2.40 ±0.04) x 10""^, in good agreement with other indicators. Error bars 
are statistical and a Salpeter IMF is assumed throughout. We find no evolution of the mass 
function in the redshift range 0.05 < z < 0.34 indicating that almost all stars were already 
formed at z ^ 0.34 with little or no star formation activity since then and that the evolution 
seen in the luminosity function must be largely due to stellar fading. 

Key words: galaxies: luminosity function, mass function - galaxies: evolution - galaxies: 
formation parameters - galaxies: statistics - galaxies: stellar content - cosmology: cosmolog- 
ical parameters 
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1 INTRODUCTION 

Recent advances in galaxy modelling codes mean that we are on 
the verge of being able to predict reliably the mass function of stel- 
lar component of galaxies in the Universe. To match this modelling 
we require a good observational determination from a large sample 
of galaxies. This has been attempted in the past by combining the 
Two Micron All Sky Surve y (2MASS) with t he 2dF Galaxy Red- 
shift Survey (2dF GRS) bv lColeetalJ <200lh and the Sl oan Dig- 
ital Sk y Survey (SDSS) Early Data Release (EDR) by Bell etjj 
i2003h . Both these methods used the K-band magnitudes of galax- 
ies in the 2MASS survey to estimate the mass of their respective 
redshift survey galaxies. 

Unfortunately neither of the redshift surveys can be used to 
their full depth when combined with the 2MASS survey - it does 
not go deep enough . The Sloan Dig i tal S ky Survey (SDSS) Dat a 
Release One (DRl. lAbazaiian et alj 1200 2). Straus s et alj J2002h ') 
gives a larger sample of galaxies (of the order 100,000 spectra) and 
it would be advantageous to be able to form a mass function from 
all of these galaxies. Obviously, just using the optical magnitudes 
from these surveys is insufficient to yield stellar mass, as some ac- 
count must be taken of AGN component, dust component, hot gas 
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and the evolution of stellar population. It would be highly unsatis- 
factory to estimate a blanket correction for these components, and 
hence an independent determination of the stellar mass is required 
for each individual galaxy in the survey. 

With the advent of efficient spectral fitting algorithms it is pos- 
sible to extract and model the stellar populations of galaxies, and 
reliably exclude the contribution to the spectrum of emission lines, 
dust and AGN. The stellar component of the mass of the galaxy 
can then be calculated, and when applied to a large enough sur- 
vey of galaxies be used to calculate the mass function of the stellar 
component. 



2 DERIVING STELLAR MASSES WITH MOPED 



2.1 MOPED 



The MOPED technique JHeavens et alj iJOOol) . iReichardt et alJ 
j200lh . |Panter et alj j2003h . lHeavens et alji2004l) ) was used to ex- 
tract the Star Formation History (SFH) of the galaxies contained 
in the SDSS Main Galaxy catalogue ( Shen et al. ( 2003)). The tech- 
nique is examined in detail in previous papers, and is equivalent to 
fitting synthetic stellar population spectra to each galaxy spectrum 
using a novel data compression and analysis algorithm. The main 
emission-line regions are excluded. Masses were calculated for a 
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Figure 1. Mass function for the SDSS DRl, with a Schechter Function fit overplotted (solid green Une). Also shown are the mass functions recovered by Cole 
et al. (2001) (dashed blue line) and Bell et al. (2003) (solid cyan line). The mass offset of the Bell data is due to a different choice of IMF. The first three bins 
have been excluded from the fit. 



complete sample inside a redshift range of 0.005 < z < 0.34 
with R band apparent m agnitudes 15.0 ^ niR ^ 17.77. Fol- 
lowing the suggestion of Ishen et alJ i2003^ we have also set a 
surface brightness limit of fiR < 23.0. Throughout we assume 
the concordant WMAP cosmology, Q,jn = 0.27, ily = 0.73, 

= 0.71 JSnergeletalJJ2003h ). 



Ho = lOO/i km s"^Mpc" 



2.2 Mass extraction 



MOPED extracts the star formation and metalicity history of a 
galaxy along with its dust content. The star formation history 
of each galaxy is modelled by 23 numbers: the star formation 
fraction in each of 11 time periods, largely spaced equally in 
log(lookback time); 11 associated metallicities of the star-forming 

[ as; a simple d ust screen characterised by an LMC extinction curve 
Gordon et alJ 12003) ) and a single dust parameter. By combining 
this information with the input stellar model spectra we can cre- 
ate a synthetic spectrum of the stellar component of the galaxy 
spectrum. MOPED adjusts the 23 parameters until the best match 
with 23 numbers derived from the galaxy spectrum is found. These 
parameters, combined with an overall normalisation of the galaxy 
spectrum, allo w the integrated ma ss created in stars to be estimated. 

Following ICole et alj i200lh to account for the mass of stars 
lost through winds and supemovae we use the recycling fraction 
R = 0.28 derived from stellar evolution theory. The final stellar 
mass of the galaxy is then ( 1 — i?) times the integrated mass of the 
different fitted populations totalled over the bins. 

The masses recovered by this method are dependent on choice 
of IMF. We have chosen a Salpeter IMF to allow direct comparison 
with theory and other predictions. Since we are using a Salpeter 
IMF with a low mass cut off at 0.1 Mq we are insensitive to any 
mass included in brown dwarfs. 

To turn the measurements of mass to a mass density we weight 



each galaxy by 1/Vmax, where Vmax is the maximum volume in 
which a given galaxy would still be included within the limits of the 
survey. In the past these have been calculated by assuming some 
luminosity evolution and calculating a K-correction from reference 
spectra. Since we have knowledge of the star formation history of 
every galaxy in the SDSS we can evolve the luminosity and surface 
brightness history of each galaxy rather than assume a fixed evolu- 
tion. To evolve the surface brightness of the galaxy we assume that 
over the redshift range studied there has been no size change of the 
galaxy. The mass estimator should be unbiased; in principle it is 
sensitive to clustering, but the effective volumes probed are large 
enough that this should be a small effect. 



To account for the three-arcsec fibre diameter we scale our 
masses up by the ratio of the flux from the photometric R band pet- 
rosian magnitudes to that from the spectroscopically determined fi- 
bre magnitudes. This is unlikely to succeed for individual galaxies, 
but there is evidence from the Petrosian and fibre colours that the re 
is no systematic offset caused overall JGlazebrook et alj y003)). 
We also find that, for galaxies of about the same mass, the interme- 
diate results such as star formation fraction show no signs of a trend 
with galaxy redshift, which one would expect if there were a signif- 
icant systematic error arising from aperture effects iJimenez et alJ 
i2004 ). 



To calculate the statistical errors on our mass function we ap- 
plied a bootstrap error algorithm. The binning operation was per- 
formed with 1000 randomly-selected resamples of the original data 
set, and the standard deviations of the heights of each peak recov- 
ered. 
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Figure 2. Mass Function with Modified Fit. Here we present a well-fitting modification to the Schechter function which gives excellent agreement with our 
results at high mass. Again, the first three bins have been excluded from the fit. The Schechter function alone is a good fit up to 10^^'''h^^MQ 



3 RESULTS 

3.1 Galaxy Stellar Mass Function 

The stellar mass function is shown in Fig.Q Between about 10® and 
10^^ h~^ Mq we find excellent agreement with results obtained by 
pre vious studies of SP SS and 2dFGRS galaxies ( Bell et al. ( 2003) 
and ICole et alJ yOOl|)), where the stellar mass is estimated more 
simply from infrared data. We are able to extend the mass range 
considerably, by around a decade in mass at the upper mass end, 
and about two decades at the lower-mass end. The stellar mass 
function of SDSS galaxies is now accurately determined between 
10^"'' and 10^^ h~'^MQ, where h is the Hubble parameter in units 
of 100km s~^ Mpc~^ . 

We fit the galaxy stellar mass function with a Schechter func- 
tion <Schechteil7T97a) ) 



(f>{Ms)dMs 






dMs 



(1) 



with best-fitting parameters (f)* = (7.7 ± 0.8) x 10"^^^ Mpc"^, 
a = -1.167 ± 0.004 and M* = (7.53 ± 0.04) x lO^^h'^MQ. 
This fit is shown overplotted in Fig|2| and is a good fit up to Ms = 
lO"-^/i-^M0. 

There is evidence for an excess over the Schechter fit at the 
high-mass end (which seems to be confirmed from dynamical mea- 
surements of the mass of SDDS galaxies, Bernardi et al. private 
communication), which can be well modelled by the addition of a 
power law over the range 11.5 < log j^g(h'^ Ms /Mq) < 12.6 of 
the form 



= <^ + 



-O" 



(2) 



with 13 = -4.03 ± 0.03 and Fc = (2.1 ± 0.2) x lO"^/;,^ Mpc"^ 
as shown in Figure |2| This excess could be due to cD galaxies or 
a failure of the modelled correction to total magnitudes for these 
extremely large galaxies. 



3.2 Evolution with Redshift 



By splitting the DRl sample by redshift we can probe the evolution 
of the stellar mass function in recent times. Fig.|3|shows the stellar 
mass function for galaxies within narrow redshift ranges. Because 
of the flux limit, there is essentially a minimum mass which can be 
probed at each redshift, but this is not a sharp cutoff because the 
galaxies have a range of star formation histories so the mapping 
from stellar mass to luminosity is not one-to-one. It is apparent 
from the figure that within the limits of the survey there is very 
little, if any, evolution in the redshift range 0.05 < z < 0.34. The 
only notable deviation from this is an apparent deficiency of high- 
mass galaxies (Ms > 10^^ H'^Mq) at very low z ~ 0.05. The 
high mass results from the lowest redshift sample should be treated 
with caution. The galaxies at the high mass end of the mass function 
are generally large in their angular size. This leads to a problem 
with "shredding" by the SDSS photometric pipeline, where large 
galaxies are treated as many smaller sources. It is thought that this 
is only really a problem for z < 0.02, but for z < 0.01 as many 
as 10% of the detections could be affected (SDSS Collaboration, 
private communication). 

The lack of evolution of the mass function with redshift is in 
contrast to the significant evolution found in the luminosity func- 
tion, where the characteristic luminosity has become fainter by 
around 0.3 magnitudes since z = 0.2, and the number density of 
bright g alaxies has de c lined b y a factor of two or more JLovedavl 
( 2003): iBlanton et alJ y003t|)). The most natural explanation is 
that the stellar mass content has hardly changed, but that the galax- 
ies have just become significantly fainter; this is expected given the 
drop-off in star formation activity, and can be illustrated by Fig.|4] 
which shows the evolution of the average stellar mass with redshift, 
for galaxy populations of different luminosities. We see clearly an 
increase in the average mass with decreasing redshift. 
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Figure 3. Mass functions for different redshift ranges. Tlie agreement is generally very good where the samples overlap, indicating that there is little evolution 
over the redshift range — 0.3. There is some discrepancy at the high-mass end in the lowest redshift range, thought to be due to Sloan photometric pipehne 
shredding large galaxy images. 
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Figure 4. The evolution of the average galaxy stellar mass for galaxies with 
(from the bottom) logio(Lfl/L0) = 10.25-10.3, 10.45-10.5, 10.65- 
10.7, 10.85 — 10.9, where Lb. is the i?— band luminosity, ii'— corrected 
using kcorrect v3_lB (Blanton et al (2003a)) 



3.3 Cosmological Baryon density in stars 

Our results can be used to give a further constraint on the con- 
tribution to the density parameter from baryons in stars, Q,h*. By 
integrating the mass over the range of the mass function we de- 
duce a value of Q.i,*h = 2.40 ± 0.04 x 10"^. This value is 
in bro a d agreeme nt with results o btained previ o usly tCole et al. 
200l [) : iBell et alj i2 003): Fukugita 'etai] h99d) : lKochanek et alj 
2001 ): lGlazebrook et aljl2003i) ). Our error is a bootstrap estimate, 
and is purely statistical; systematic errors such as the choice of IMF 
have not been included. 



4 DISCUSSION AND CONCLUSIONS 

We have calculated the stellar mass function of 96,545 galaxies 
in the Data Release of the Sloan Digital Sky Survey over about 

5 decades of mass. The results are in good agreement with previ- 
ous studies, where the stellar masses were estimated more simply 
from infrared data. The range probed is considerably extended, and 
differences in derived parameters are most likely due to different 
assumed initial mass functions. 

Also of importance is the difference in numerical size and red- 
shift range of the sample used to generate the stellar mass function, 
since we are not restricted to galaxies which appear in the range of 
galaxies in 2MASS. 

In contrast to the luminosity function, the mass function shows 
no evidence for evolution with redshift. The luminosity function 
shows a fading of the characteristic luminosity, by a factor of 
about 1.35 since z = 0.2 (Lovedav (2003); Blanton et al. (2003b)). 
The simplest interpretation is that the galaxy stellar masses do not 
evolve significan tly (and this is supp orted by the star formation 
rates reported bv lHeavens et alJ yOO^), but that individual galax- 
ies fade. Fig.|4]shows the average mass of massive galaxies in small 
luminosity redshift ranges, as a function of the observed redshift. 
We see the expected trend: massive galaxies are typically brighter at 
higher redshift, scaling roughly as logjQ Ms{z) ~ logig Ms{z — 
0) — OAz. This represents a typical fading of around 20% over the 
redshift range z = 0.2 to the present. 
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